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We present extensive first principles density functional theory (DFT) calculations dedicated to 
analyze the magnetic and electronic properties of small V„ clusters (n=l,2,3,4,5,6) embedded in a 
Cu fee matrix. We consider different cluster structures such as: i) a single V impurity, ii) several 
V2 dimers having different interatomic distance and varying local atomic environment, iii) V3 and 
iv) V4 clusters for which we assume compact as well as 2- and 1-dimensional atomic configurations 
and finally, in the case of the v) V5 and vi) V6 structures we consider a square pyramid and a 
square bipyramid together with linear arrays, respectively. In all cases, the V atoms are embedded 
as substitutional impurities in the Cu network. In general, and as in the free standing case, we 
have found that the V clusters tend to form compact atomic arrays within the cooper matrix. Our 
calculated non spin-polarized density of states at the V sites shows a complex peaked structure 
around the Fermi level that strongly changes as a function of both the interatomic distance and 
local atomic environment, a result that anticipates a non trivial magnetic behavior. In fact, our 
DFT calculations reveal, in each one of our clusters systems, the existence of different magnetic 
solutions (ferromagnetic, ferrimagnetic, and antiferromagnetic) with very small energy differences 
among them, a result that could lead to the existence of complex finite-temperature magnetic 
properties. Finally, we compare our results with recent experimental measurements. 



PACS numbers: 75.75. +a, 73.22.-f, 75.50.-y 

I. INTRODUCTION 

The discovery of magnetism in V metal fine particles 1 
with sizes and temperatures ranging from 90 to 300 A 
and from 4.2 to 100 K, respectively, served as a starting 
point for a large number of theoretical and experimen- 
tal studies of magnetic properties in various types of V 
low-dimensional systems^- 2 . In particular, bimetallic V- 
based materials such as bulk alloys^, thin filmsi2, and 
deposited cluster o 10 ' 13 have attracted considerable atten- 
tion because of their expected novel optical properties, 
catalytic activity, and magnetic behavior. 

In general the interaction between two components in 
magnetic bimetallic nanostructures introduces a mutual 
influence on neighboring atoms and leads to the unique 
properties already reported for these kinds of materials. 
Furthermore, it is also clear that their measured aver- 
age magnetic behavior is not necessarily given by the 
average properties of their corresponding isolated con- 
stituents, and that the observed phenomena are expected 
to strongly depend on the precise details of the local ge- 
ometrical and chemical environment. 

As is well known, V in the solid state is nonmagnetic; 
however, it has been clearly established that it could ex- 
hibit a magnetic moment under certain conditions like 
loss of coordination, when its atomic volume is increased, 
or when it is found in the presence of 3d ferromagnetic 
materials, due to its large paramagnetic susceptibility. 
Interestingly, in the last years it has been also demon- 
strated that V can become magnetic when it is associ- 
ated with other nonmagnetic elements such as copper— 
or golcU^. The previous findings are highly relevant since, 



on the one hand, having two nonmagnetic elements that 
when combined give raise to magnetic order opens new 
possibilities to improve the magnetic properties of mate- 
rials and, on the other hand, from the theoretical point 
of view it offers a unique opportunity to analyze the role 
of the structure and dimensionality on the spontaneous 
appearance of local magnetic moments in metallic nanos- 
tructured materials. 

Huttel et al^ have found experimentally the existence 
of a magnetic moment in V atoms embedded in a Cu 
matrix. In particular, they have analyzed samples hav- 
ing a CU96.7V3.3 composition by means of x-ray absorp- 
tion spectroscopy and x-ray magnetic circular dichro- 
ism experiments. The samples were prepared by a co- 
evaporation process of both V and Cu over a Cu sub- 
strate and the as-prepared Cu-V samples were capped 
with a few monolayers of Cu in order to avoid oxidation 
of the surface alloy. Interestingly the authors especulate 
that, depending on the deposition temperature, different 
V species seems to be formed (such as isolated V im- 
purities and V particles of different sizes) and that the 
magnetic properties of the alloy strongly depend on the 
presence and concentration of these species. They also 
showed that the V magnetic moment depends on the for- 
mation temperature of the samples and found that the 
alloys prepared at lower temperatures are more magnetic. 
However, despite all the previous extensive experimen- 
tal characterization and reported trends, there is still no 
clear information concerning the precise micro-structural 
features present in the alloys, neither a clear explanation 
for the temperature dependence of the average magneti- 
zation observed in the samples. 



From the point of view of theory, it is important to 
mention that the possible existence of structural changes, 
disordered growth, and clustering of the V species in the 
samples makes more difficult the theoretical analysis of 
the measured data. However, we believe that the Cu- 

V alloys synthesized by Huttel and co-workers can be 
approximately modeled by an infinite Cu matrix with 

V atoms and clusters embedded in it. In addition, we 
expect that the influence of the Cu matrix on the V 
atoms will be mainly due to the hybridization of the V 
(i-electrons with the s-levels of the neighboring Cu atoms 
and, in this sense, a local analysis in which the behavior 
of the spin moments at the V sites is studied for different 
geometrical and chemical environments should be very 
important and needs to be performed in order to shed 
some more light in to the possible origin of the measured 
magnetization data. 

In this work, we present a systematic first principles 
theoretical study of the magnetic properties of model V- 
dilute fee Cu-V alloys by using the density functional 
theory (DFT) approach together with the pseudopoten- 
tial approximation for the electron-ion interaction. All 
the calculations are performed with the Plane Wave Self 
Consistent Field (PWscf) codei£. We put special em- 
phasis on: i) the electronic and magnetic properties of 
isolated atoms as well as of small clusters embedded in a 
fee Cu matrix, and ii) the possible existence of multiple 
magnetic solutions in the systems, since if small energy 
differences are found among them magnetic fluctuations 
could be present in the samples. 

The remainder of the paper is organized as follows. 
In the next section, the theoretical model used for the 
calculations is briefly explained. In Sec. IIIII results for 
the atomic spin magnetic moments, average spin mag- 
netic moments, and differences in total energies (DTEs) 
for a variety of V^v clusters embedded in a Cu fee ma- 
trix with N < 4 atoms, and some special cases of TV = 5 
(square pyramide) and N = 6 (square bipyramide), as 
well as infinite linear chain are presented and discussed. 
For each system we have studied the nonmagnetic (NM), 
ferromagnetic (FM), and antiferromagnetic (AF) solu- 
tions, i.e, we have used these magnetic configurations as 
inputs in our calculations. Finally, Sec. IIVI summarizes 
the main conclusions. 



II. 



THEORETICAL MODEL 



The PWscfi£ code is based on the DFT. The pseu- 
dopotentials used here were generated using Vanderbilt 
codeii with scalar relativistic calculations and non linear 
core corrections. The calculations, which were restricted 
to fixed atomic positions (relaxation was not allowed), 
were performed with the generalized gradient approxima- 
tion of Perdew-Burke-Ernzerhof (GGA-PBE) 18 exchange 
correlation functional. The calculations were performed 
using a fee structure with Cu lattice parameter (3.61 A 
) and the cell containing 48 (54) atoms (see Fig. Q]), or- 



(a) Top view for the eell with 48 atoms 
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(c) Side view for cell with 48 atoms 



(b) Top view for the cell with 54 atoms 
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(d) Side view for cell with 54 atoms 



FIG. 1. Schema of the crystal cell used, (a) Top view for the 
crystal cell with 48 atoms, (b) top view for the cell with 54 
atoms, (c) and (d) are side view for the crystal cell with 48 
and 54 atoms, respectively. The dotted line is the crystal cell 
and the a is the Cu lattice parameter (3.61 A ). The empty 
circles represent the Cu or V atoms. 



dered as follows: six planes with 8(9) atoms per plane for 
the impurity, dimer, trimer and tetramer (pyramide and 
bipyramide), and the V atoms are located at the ideal 
positions of the fee crystal. The Monkhorst-Pack scheme 
was used to define the k points and the calculations for 
each cell were performed with a grid of 4x4x4 mesh in 
Af-space. We have done some tests with a grid of 5x5x5 
and 6x6x6 mesh, and the values are only 4% less than 
those presented here. A cutoff energy of 35 Ry was used 
for the plane waves expansion of the pseudowave func- 
tions (560 Ry for the charge density and potential). 

Given that our calculations were restricted to fixed 
atomic positions, some words about the relaxation ef- 
fects in this kind of systems seem to be in order. On 
the one hand, Ramanathan et aS^, who have studied 
the magnetism of a V overlayer on Nb(001), have found 
that the magnetism is drastically reduced when the re- 
laxation effects are considered. On the other hand, We- 
ber et al 1 ^, who have studied the magnetism of V clus- 
ters supported on Cu(001), have analysed the effect of 
relaxation qualitatively, and they have concluded that 
relaxation does not yield any significant change in the 
magnetic moments. Following that, we performed cal- 
culations for the V impurity and dimer using a cell 
with 54 atoms, first without considering relaxation and 
then considering partial relaxation (only nearest neigh- 
bors could reach the equilibrium positions) and we have 
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FIG. 2. Illustration of the lattice structure of Vjv clusters embedded as substitutional impurities in bulk Cu (fee structure). 
The cluster atoms are represented by filled circles and the Cu host atoms by open circles. The numbers label the nonequivalent 
atomic sites i. — (ifm/^af refers to the magnetic moments (in fis) for the ferromagnetic (FM) and antiferromag- 

netic (AF) solutions, respectively. DTE=DTEjva//DTEfa//DTEaf refers to the differences in total energies (DTEs) (in 
meV/Atom) for the nonmagnetic (NM), FM, and AF solutions, respectively. In the special case of the square bipyramide 
[Fig. [2] (y)] we include two kinds of antiferromagnetic solutions, namely, AF1 and AF2, whose initial magnetic configurations 
are, respectively, (1 t, 2 I, 3 t, 4 I, 5 f, 6 t) and (I f, 2 |, 3 f, 4 |, 5 f, 6 |). In this case the DTEs are presented as follows: 
DTE=DTEjvm/DTEfa//DTEafi/DTEaf2- The ground state energies (denoted by 0.00) refer to the minimum energies out 
of Vjv clusters with the same TV. 



obtained that the differences in percentage terms be- 
tween non-relaxation and partial relaxation calculations 
are AE < 4% ; Ad =~ 1.0% and A/i =~ 2% for the 
difference of total energies, difference of distances, and 
difference of magnetic moments, respectively. 

III. RESULTS AND DISCUSSION 
A. Geometrical construction 

The V clusters geometric arrangements (as well as the 
lattice structure around them) that we have studied are 
schematically illustrated in Fig. [2J The V embedded 



clusters are located as substitutional impurities in the Cu 
fee matrix. We present the atomic magnetic moments as 
n(i) = ^fm/^af in Hb for the FM and AF solutions, re- 
spectively. We also present the difference in total energies 
as DTE = DTE N m/DTE F m/DTEaf (in meV/atom) 
for the NM, FM, and AF solutions, respectively. In 
the special case of the square bipyramide [Fig. [5] (y)] 
we present the difference in total energies as following 
DTE = DTE NM /DTE FM /DTE AF1 /DTE AF2 / where 
AF1 and AF2 stand, respectively, for the next initial 
magnetic configurations: (1 f, 2 J,, 3 t, 4 J,, 5 J) 6 f) and 
(1 t> 2 |,3 t) 4 1,5 f, 6 i). The ground states energies 
(denoted by 0.00) refer to the minimum energies out of 
Vjv clusters with fixed N. In all the cases below, a = 3.61 
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A is the Cu lattice parameter. Fig. [2] (a) represents the 
V impurity. The V dimers are also in the first row: (b), 
(c), (d), (e), and (f). The respective interatomic dis- 
tances are: a /y/2, «^3/2, a, ay/3, and ay/2. The V 
trimers are in the second row. (g) represents the small- 
est equilateral triangle whose sides measure a/y/2; (h) 
an isosceles right angle triangle whose hypotenuse and 
legs measure, respectively, a and a/y/2; (i) a right angle 
triangle whose hypotenuse measures ay^3/2, and whose 
legs measures a/y/2 and a; (j) a linear trimer whose in- 
teratomic distances are equal to a/y/2; (k) an isosceles 
right angle triangle whose hypotenuse and legs measure, 
respectively, ay/2 and a; (1) a right angle triangle whose 
hypotenuse measures ay/3, and whose legs measure a and 
ay/2; and (m) an equilateral triangle whose sides measure 
ay/2. The V tetramers are in the third and fourth rows, 
(n) represents a regular tetrahedron whose interatomic 
distances are a/y/2; (o) a regular rhombus whose sides 
measure a/y/2; (p) a square whose sides measure a/y/2; 
(q) a triangular tetramer with two bonds of length equal 
to a/y/2 and two bonds with length equal to a; (r) a par- 
allelogram whose sides measure ay/3/2 and a/y/2; (s) a 
continuous linear tetramer whose interatomic distances 
are equal to a/y/2; (t) a triangular tetramer with two 
bonds of length equal to a/ y/2, one bond of length equal 
to a, and one bond of length equal to ay/3; (u) a rectan- 
gle with sides equal to a/y/2 and a; (v) a square with side 
equal to a; (w) a rectangle with sides equal to ay/2 and 
a; the special cases are (x) a square pyramid whose base 
measures a/y/2 x a/y/2 and whose sides are all equilat- 
eral triangles of edge-length parameter equal to a/y/2; (y) 
a square bipyramid built from two pyramids connected 
base to base; and finally (z) a continuous linear chain 
whose interatomic distances are equal to a/ \/2 

Fig. [2] displays the well-known richness and complex- 
ity in the magnetic behaviour characterizing the 3d tran- 
sition metal low-dimensional systems. There we can 
see four kinds of magnetic orders: ferromagnetic, fer- 
rimagnetic, antiferromagnetic, and nonmagnetic, which 
depend in a very complicated way on the interatomic 
distances, the geometrical structures, and the number of 
atoms in the clusters. 



B. Electronic behavior 

Let us first briefly consider the electronic properties 
derived of our calculations on Vn clusters embedded in 
a Cu fee matrix. Fig. [3] illustrates the d-component of 
the nonmagnetic local densities of states (NMLDOS) for 
all the V dimers and the tetrahedron. The blue, red, 
violet, black, and green lines represent, respectively, the 
NMLDOS for the V dimers illustrated in Fig. [1(b), (c), 
(d), (e), and (f); while the orange line represents the 
NMD OS for the tetrahedron [Fig. [2] (n)]. To study the 
electronic properties of these systems we use the Stoner 




E - EF (eV) 

FIG. 3. The nonmagnetic local density of states (NMLDOS) 
for V atoms (only d component). The blue, red, violet, black, 
and green lines represent, respectively, the NLDOS for the 
dimers illustrated in Fig. [2](b), (c), (d), (e), (f). The orange 
line represents the NLDOS for the tetrahedron tetramer [Fig. 
[2](n)]. The dashed line represents the Fermi level. The inset 
shows a magnified view of the LDOS near the Fermi level. 

criterion. This is given by the following relation: 

S = I S N NM (E F ) 

being S the Stoner factor, 1$ the Stoner parameter and 
N NM (E F ) the value of the NMLDOS at the Fermi level. 
The I s for V is equal to 0.3536 eV-22. The Stoner crite- 
rion states that the nonmagnetic state will be unstable 
to magnetism if S > 1, while the nonmagnetic state will 
be more favourable if S < 1 . 

The behaviour of the NMLDOS in Fig. [3] is similar to 
that reported by Liu, Khana, and Jena^i, who studied 
the magnetism of free small V clusters, in the sense that 
near the Fermi energy the NMLDOS splits into two peaks 
with the Fermi energy centered in the valley. This split- 
ting is remarkably pronounced in the case of the tetrahe- 
dron. We can see that for all the V dimers the value of 
the NMLDOS at the Fermi level are high enough (greater 
than I /I s = 2.8280) to get a Stoner factor S > 1; while 
for the tetrahedron tetramer the NMLDOS at the Fermi 
level is not (it is much less than 1/Is = 2.8280). So from 
the Stoner criterion alone we can predict magnetism for 
all the V dimer and absence of magnetism for the tetra- 
hedron tetramer. The greatest value of the NMLDOS at 
Fermi level corresponds to the dimer whose interatomic 
distance is ay/2 — 5.1 A (see [2] (f) and green line in 
Fig. which, according to the values of the magnetic 
moments shown in Fig. [2] is that with the greatest mag- 
netic moment out of the V dimers. Additionally, from 
the NMLDOS we can anticipate a nonlinear behaviour 
of the magnetic moments of V dimers as a function of 
the interatomic distances. 
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C. Magnetic properties 

Let us now consider the magnetic behaviour of each 
cluster as a function of its size. For the V impurity, the 
magnetic moment is 1.87 /j,b, which is less than that of 
the V adsorbated atom on Cu(OOl) surface reported by 
Weber et (3.03 /is) and greater than that of the bulk 
impurity reported by Bliigel et aS^ (1.10 \xb)- 

For the V dimers, we can see in Fig. H](b)-(f) that the 
magnetic states are energetically ordered in such a way 
that E NM > Efm > Eaf (Enm > Eaf > E F m) pro- 
vided that the interatomic distance d < a (d > a). So, in 
the first (second) case the V atoms that form the dimers 
couple antiferromagnetically (ferromagnetically) . In any 
case, for the V dimers, it is energetically easier to show 
some kind of magnetic coupling than to become nonmag- 
netic. Interestingly, the V dimer ground state does not 
correspond to the V dimer with the minimum interatomic 
distance, namely a/y/2, but to that with the third short- 
est one: a\Ji/2 [see Fig. [2] (b) and (c)]. The magnetic 
moment per V atom for the ground state V dimer is equal 
to ± 1.62 fig which is less than adsorbate dimer mag- 
netic moment reported by Weber et al~2- (2.85 hb), but 
greater than the bulk dimer magnetic moment reported 
by Bliigel et al^ (0.00 ^s). The first three excited states 
correspond to the V dimer with d = aj \[2. The DTEs for 
these states are equal to 0.8, 1.0, and 1.51 meV/atom, 
respectively, for the antiferromagnetic (/i = ±1.27/is), 
ferromagnetic — 1.00/is), and nonmagnetic solutions. 
So, from the ground state, the most probable transition 
is that in which the V dimer reduces both its interatomic 
distance and the magnitude of its magnetic moment, and 
keeps the same kind of magnetic coupling (antiferromag- 
netic in this case). Note that it is more probable that 
from the ground state the V dimer reduces its interatomic 
distance than it changes its coupling from antiferromag- 
netic to ferromagnetic, since for the first (second) transi- 
tion it needs to increase its energy in 0.8 meV/atom (1.6 
meV/atom). For the minimum interatomic distance, a 
transition from antiferromagnetic to ferromagnetic cou- 
pling can occur if the V dimer increases its energy in 
0.2 meV j atom. This transition implies a reduction of 
the magnitude of the magnetic moment from 1.27 \lb to 
1.00 fiB- In general, the magnitudes of the magnetic mo- 
ments for the ferromagnetic solutions are less than those 
for the antiferromagntic ones (the only exception is that 
of the V dimer with d = ay/2 in Fig. [2] (f)). Also, in 
general, the magnitude of the magnetic moments tend to 
increase (although not monotonously) as the interatomic 
distance increases. However, in Fig. [4] we can observe 
that the magnitude of the magnetic moments increase as 
the DTEs increase in a strictly monotonous way. Note 
also that although the interatomic distance for the V 
dimer represented in Fig. [2] (e) is greater than that for 
the V dimer represented in Fig. [2] (f) its magnetic mo- 
ment is smaller (1.85 fis vs. 2 fis, and 1.91 [Ib vs 1.95 
/iB, for the FM and AF solutions, respectively) maybe 
due to the presence of a Cu atom between the V atoms. 



Finally, to compare our calculations with the measure- 
ments made by Huttel et al^, we note first that for the 
V dimer the concentration [Cu 4 6/48V 2 /48 (Cuo.96Vo.04)] 
coincides approximately with that reported by Huttel et 
al*£ (Cuo.967Vo.033) who reported a mean magnetic mo- 
ment in V atoms of approximately 1.00 /is, which nicely 
agrees with the V dimer first state (FM). 

In general, the trimers couple ferromagnetically (the 
only exceptions are the linear chain in Fig. [2] (j) and 
the isosceles triangle in Fig. [2] (k)). The order of the 
magnetic states is not so clearly defined in terms of the 
geometric characteristics of the clusters as it was in the 
case of the V dimers. As in the case of the V dimers, 
for the V trimers, it is energetically easier to show some 
kind of magnetic coupling than to become nonmagnetic. 
The V trimer ground state corresponds to the more com- 
pact equilateral triangle in the Fig. [2] (g), whose inter- 
atomic distances are equal to a/y/2. Practically, this 
ground state has 2-fold degeneracy since the difference 
in the DTE for two states of lowest energy is only 0.05 
meV/ atom. So, we can say that there the ferromagnetic 
and nonmagnetic states coexist. For the ferromagnetic 
solution, the magnetic moment per V atom is equal to 

0. 46 fiB- If we think of this trimer as a superior stage in 
the evolution of the dimer in Fig. [5] (b) as the number 
of atoms increase, we note, as expected, a diminishing in 
the magnetic moment as a consequence of the increment 
in the coordination number. For the V trimers in Figs. 
El (g)-(j)i it is easier to change their magnetic coupling 
than to change their geometric configuration. The op- 
posite occurs for the V trimers in Figs. [5] (k)-(m), for 
example, if the V trimer in Fig. [5] (k) is in its ferromag- 
netic state, it costs less energy to change its geometric 
configuration to that shown in Fig. [2] (1) than to become 
nonmagnetic while keeping the initial geometric configu- 
ration. Finally, we note that for the linear trimer (Fig. 
[2](i)) the magnetic moments are 1.51 /is, —1.08 /.ib and 
1.51 hb, these are lower than the moments found by We- 
ber et al*£ in the adchain of 3 V atoms on Cu(001) fee. 

It is not easy for the V tetramers either to give general 
rules about their magnetic behaviour in terms of the ener- 
gies of their states. The ground state corresponds to the 
triangular pyramide. This state has 3-fold denegeracy, 

1. e., there coexists the nonmagnetic, ferromagnetic, and 
antiferromagnetic states. In any case, the magnetic mo- 
ments per V atom is equal to zero [see Fig. [2] (n)]. This 
does not mean, however, that the V tetramer represents 
a critical size of V clusters from which the magnetism 
vanishes since as it will be seen below, V clusters with 
five and six atoms give nonzero magnetic moments. So, 
according to our calculations the V magnetic moments 
fall to zero more slowly than those reported by Bliigel 
et al.—. In general, we note that for small V clusters 
embedded in a Cu fee matrix the most stable structures 
trend to form compact arrangements: the smallest equi- 
lateral trimer and the triangular pyramide. The excep- 
tion is the V dimer, whose ground state occurs for the 
dimer with the third shortest interatomic distance. For 
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FIG. 4. Mean magnetic moments as a function of the dif- 
ference of the total energy. The black, red and blue lines are 
related to V dimers, trimers, and tetramers, respectively. The 
dashed lines represent the corresponding linear regressions. 

the tetramers, the first excited state corresponds to the 

V rhombus in Fig. [2] (o). This state corresponds to the 
ferromagnetic solution but in reality the V atoms cou- 
ple ferrimagnetically and has a mean magnetic moment 
of 0.31 hb- The linear chain [Fig. [5] (s)] shows a ferri- 
magnetic order. As expected, its magnetic moments are 
smaller than those corresponding to the V dimer and the 

V linear chain trimer. 

Beyond TV = 4, we have considered three other sys- 
tems, namely, square pyramide (TV — 5) and square 
bipyramide (TV — 6), and continuous linear chain. The 
ground state for the square pyramide corresponds to the 
ferromagnetic solution with a mean magnetic moment 
per V atom equal to 0.70/is. For the square bipyra- 
mide we have considered four solutions, namely, non- 
magnetic, ferromagnetic, antiferromagnetic one (AF1), 
and antiferromagnetic two (AF2). The input for AF1 
and AF2 are, respectively: (1 f, 2 4, 3 f, 4 4, 5 f, 6 f) and 
(1 ti 2 |,3 fj 4 4,5 fj 6 4) where the numbers label the 

V atomic sites and the arrows the directions of the input 
spin magnetic moments. The ground state corresponds 
to the ferromagnetic solution with a mean magnetic mo- 
ment per V atom equal to 0.39 /j,b- Practically, the AF1 
and AF2 solutions are degenerate. For continuous linear 
chain, we have studied both ferromagnetic and antifer- 
romagnetic solutions. The ground state is the antifer- 
romagnetic solution with magnetic moments at V sites 
equal to ±1.16/xb. If consider the evolution of the linear 
chain from the dimer to the infinite linear chain keeping 
the same interatomic distance a/V2 [see Fig. [2](b), (j), 
(s), and (z)], we note the V atomos always couple anti- 
ferromagnetically. Besides, the V atoms located at the 
extreme sites of the chain trend to the limit of the infi- 
nite chain more slowly (1.27, 1.51, 1.39, 1.16 /is) than 



the V atoms located at inner sites of the chain (1.08, 1.17, 
1.16 (J-b)- We think that for the square pyramide, square 
bipyramide, and infinite linear chain an increase of the 
degeneracy can be obtained, because it is possible to find 
other magnetic states (antiferromagnetic by blocks, for 
example). For TV > 7, we also think that it is possible to 
find magnetic degenerate states with magnetic moments 
per V atom eventually diminishing to zero (for TV large 
enough) . 

In the Fig. @]we illustrate the mean magnetic moments 
as a function as the DTE, only for the FM state. For 

V dimers, the mean magnetic moments increase as the 
DTE increases. We can say that for room temperatures, 
it will be possible to find higher magnetic moments. For 

V trimers and tetramers, the behaviour is more compli- 
cated, even so we can observe that the mean values show 
a positive slope. 

IV. CONCLUSIONS 

We conclude our work as follows: i) we have studied 
the magnetic behaviour of small V clusters (TV = 1,2,3 
and 4 atoms) and some other special geometries for 
TV = 5 (square pyramide), and TV = 6 (square bipyra- 
mide), as well as an infinite linear chain embedded in a 
Cu fee matrix using an ab initio code; ii) we have repro- 
duced theoretically an experimental result by Hutted for 
a V concentration of approximately Cuo.96Vo.04; Hi) for 
the ferromagnetic solutions, the magnetic moments per 

V atom decrease as the cluster size increases; this mag- 
netic behaviour is in agreement with previous theoretical 
studie o 13 ' 23 ; iv) the magnetic moments per V atom have 
a strong dependence of the geometric structure^ 3 -; v) the 
mean magnetic moment increases with the difference of 
total energy; vi) compact structures have lower energy 
(they have lower mean magnetic moment too), an ex- 
ception being the V dimer which prefers the third inter- 
atomic distance; whereas extended structures have both 
higher energy and mean magnetic moments; vii) for each 
one of our considerd cluster systems we have found the 
existence of different magnetic solutions (nonmagnetic, 
ferromagnetic, and antiferromagnetic) with very small 
energy differences among them, a result that we believe 
could lead to the existence of complex finite-temperature 
magnetic properties. 
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